Mechanism of formation of half-doped stripes in underdoped cuprates 
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Using a variational Monte-Carlo approach with a recently proposed stripe wave function, we 
showed that the strong correlation included in a t — J— type model has essentially all the necessary 
ingredients to form these stripes with modulations of charge density, spin magnetization, and pair 
field. If a perturbative effect of electron-phonon coupling to renormalize the effective mass or the 
hopping rate of holes is considered with the model, we find the half-doped stripes, which has on 
the average one half of a hole in one period of charge modulation, to be most stable, energetic wise 
in the underdoped region, 1/12 < 5 < 1/8. This is in good agreement with the observation in the 
neutron scattering experiments. We also find long range Coulomb interaction to be less effective in 
the formation of half-doped stripes. 

PACS numbers: 71.10.Fd, 71.38.-k, 75.50.Ec 



The mechanism for the formation of the stripe states 
has been one of the most important issues in understand- 
ing high-temperature superconducting cuprates. Al- 
though many theoretical works based on the electronic 
model alone have succeeded in comprehending experi- 
mental results qualitatively [l], Q , there are many puz- 
zles on the stripe formation in the underdoped cuprates 
La 2 - x Sr x Cu0 4 (LSCO) and La 2 ^Ba x Cu0 4 (LBCO) 
with 1/8 doping most stable @, |J H, lE 0)- For in- 
stance, the doping dependence of the incommensurate 
magnetic peaks associated with the stripe determined 
by neutron scattering experiment on LSCO obeys the 
so-called "Yamada plot" 3] where it points to the exis- 
tence of the half-doped stripe with average of 1/2 hole 
in one charge modulation period below 1/8 hole den- 
sity. Furthermore, at 1/8 doping of LBCO, the local 
density of states still has a V-shape and a node at low 
energy even though the superconductivity is almost 
completely suppressed by the static stripe. Surprisingly, 
a recent scanning tunneling spectroscopy has observed 
cluster glasses with randomly oriented domains of stripe- 
like patterns in two other families of Ca 2 - X Na x Cu02Cl 2 
and Bi2Sr2Dyo.2Cao.sCu208+s Q- So far all these re- 
sults are yet to be explained theoretically in a strongly- 
correlated framework. 

Many early theoretical works have found that stripe- 
ordered formation may be present in the extended t — J 



Hamiltonian [l(| U, 12, [l3|. In particular, using the 
density matrix renormalization group (DMRG) in the 
standard t — J models, the mechanism of the half-doped 
stripe appears to have been understood by antifcrromag- 
nctic correlation across the hole (hfj or optimum energy 
per hole of a domain wall [l5| . However due to small lat- 
tice size, weak coupling treatment, and special boundary 
conditions used in these works, this issue is still being 
debated [lj|. The result of DMRG indicated the stripe 
to be most stable when the second nearest-neighbor hop- 
ping, t', is very small. This result is exactly opposite 
to the previous variational Monte Carlo (VMC) calcu- 
lation 12j. Based upon the idea of d-w&ve resonating- 



valence-bond (<i-RVB) state, recently we have proposed a 
new stripe-ordered wave function which involves modula- 
tions of charge, spin magnetization, and pair field. These 
stripe states either in a periodic pattern or a randomly 
oriented domain have almost the same energy as the uni- 
form d-RVB state for almost all values of t'/t in the ex- 
tended t — J model 17J. This unexpected energy degener- 



acy of stripe-ordered states seems to be quite compatible 
with the proposal of fluctuating stripes [18j . However, 
this degeneracy between different stripe orders and a uni- 
form state in the t — J— type model fails to explain the 
prominence of the half-doped stripes. 

There are several possibilities for the formation of 
stripes. Long-range Coulomb interaction between holes 
has been proposed sometime ago to transform the frus- 
trated phase separation into a stripe state [SHI. On 
the other hand, a number of experiments have presented 
evidences that phonon or lattice effects are apparently 
present in these materials [U [22], [2^, 24 1 . Particularly, 



due to the structural instability, the isotope effect in 
LSCO and LBCO is also peculiarly strong near 1 /8 dop- 
ing [25||. Moreover, to understand a kink structure in 
angle-resolved photoemission spectrum [26[, many the- 



oretical works have accepted that the electron-phonon 
cou plings s hould be considered in the t — J— type models 

In this Letter, we shall consider both long-range 
Coulomb interaction and electron-phonon interaction on 
the formation of stripes, in particular half-doped stripes, 
within the context of strong correlation as in the ex- 
tended t — J model. However electron-phonon interac- 
tion has many effects. We will only examine the simplest 
effect of mass renormalization of charges due to phonon 
coupling when the coupling strength is rather weak. The 
renormalization effect depending on the local charge den- 
sity is treated self-consistently in a VMC method which 
takes into account the strong correlation accurately. In 
the latter half of the paper, we present results of includ- 
ing long-ranged repulsive interaction between holes in the 
t — J Hamiltonian without electron-phonon interaction. 
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Let us begin with an extended t — J Hamiltonian in a 
two-dimensional square lattice which is defined as, 



H 



\ ZCT 



<ij> 



s 3 . 



(1) 



The hopping amplitude 



i, i', and f" for sites 



and j being the nearest-, the second-nearest, and the 
third-nearest-neighbors, respectively. Other notations 
are standard. In the following, we focus on the case 
(t',t",J)/t = (-0.2,0.1,0.3) as we have shown [HI that 
comparing with the uniform state, energies of the stripe 
state are not very sensitive to the values of t'/J and 
t" / J. Following our treatment in Ref. 17| , the stripe wave 
function is constructed by the mean-field Bogoliubov de 
Gennes method. We assume that the spatial dependence 
of the local charge density pi, staggered magnetization 
mi, and nearest-neighbor pair field Ay has the simple 
forms: 



Pi = p v coa[2q ■ (j/, - y )], 



n%i = m„sin[<? • (yi - y )], 



(2) 



(3) 



A 
A 



= Af cos[2g • ( Vi - y )} 

, M 



A c , 



Af cos[2g • (yi - y ) +q] + A%, (4) 



where q = ir/a c and a c is the period of charge modu- 
lation. Here the stripe extends uniformly along x di- 
rection, and yo is 1/2 (0) for the bond-centered (site- 
centered) stripe. Sited-centered stripe will not be consid- 
ered throughout the paper unless specifically mentioned. 
This state called AF-RVB stripe state 17| is quite dif- 
ferent from the antiphase stripe studied by others [12]. 
The magnitude of pair field and the magnetization are 
larger at sites with smaller hole density. Once the vari- 
ational parameters are given, we diagonalize the mean- 
field Hamiltonian. Then, the trial wave function with a 
Gutzwiller projector P can be constructed by creating 
all negative energy states 7^ and annihilating all positive 
energy states y n on a vacuum of electrons |0), 



PPjPn, 



n 



(5) 



with the fixed number of electrons N„. Here we introduce 



a particle-hole transformation [32j in Eq. ([5]) , 



J 



fi and 



d\. Using this transformation, we avoid the pos- 
sibility of having a divergent determinant because of the 
presence of nodes in RVB-type wave function with the 
fixed number of particles using periodic boundary con- 
dition. We also introduce a hole-hole repulsive Jastrow 
factors Pj in order to greatly imp rove the variational en- 
ergies (See the details in Ref.|l7j]). 

The simplest effect of an electron-phonon interaction 
is to renormalize the mass of charge carriers. To take 
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FIG. 1: (Color online) The optimized energy difference be- 
tween the AF-RVB stripe and uniform d-RVB states vs the 
doping density S for different periods of stripe. 10ao is calcu- 
lated in 20 x 20 lattice and others in 24 x 24 one. Here A is 
0.25. Inset: the electron-phonon coupling constant vs hole 
density. Filled circles show the data points from Ref.[3l|. The 
dashed line is 1 — 35. 



into account of it within the extended t — J model, 
we assume the hopping terms £y in Eq.fTJ, which are 
inversely proportional to the effective mass, are mod- 
ified due to the spatial variation of hole density. We 
neglect the renormalization of the exchange interaction 
J since it is only dominant in the very undcrdoped re- 
gion with a very small effect [28| . We also do not con- 
sider the on-site e nerg y that may arise from the electron- 
phonon coupling [33j. If the on-site energy varies lin- 
early with the local hole density, it will have the same 
energy irrespective of the period of periodic stripes. Ac- 
cording to the experimental results that the electron- 
phonon coupling strength A is reduced by the hole den- 
sity S (3l| . we assume a linear relation between A and 5, 
A(<5)/A(0) = f(S) = 1-3(5 (See the inset in FigHJ. Ac- 
cordingly, sites with a larger hole density in a modulated 
charge pattern have larger iy [34| . Thus Uj in Eq. (p} is 



renormalized to given by 



1 - A 



(6) 



where n\ is the hole density at site i and A = \(0)/ir. 
We estimate a reasonable range of A should be from 
to 0.25 Since A depends on the hole density n l h , we 
have to do a self-consistent calculation to find the mini- 
mal variational energy. Here we use the iterative method 
to achieve the convergence. First we assume an initial set 
of values for the hole density at different sites such that 
t*j are given, then we can proceed to optimize the trial 
wave functions. If the hole densities calculated by this 
optimized wave function do not agree with the initial in- 
put values, we will then use the calculated density, hence 
its t*j , as input value for the next round of optimization. 
This process is repeated until we have a converged result. 
Sometimes we also try to use either the overestimated or 
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FIG. 2: (Color online) (a) The optimized energy difference 
between the AF-RVB stripe and uniform d-RVB states vs 
A. (b) The average hole density of the AF-RVB stripes with 
different charge periods vs A. For 6ao and 8ao, empty [filled] 
symbols indicate rih(y = 1) [nh(y = 2)]. As for 12oo, empty, 
half-filled, and filled symbols correspond to nu(y = 1), rih{y = 
2), and nh(y = 3), respectively. Here the doping is 1/8. 



underestimated densities as initial input values to exam- 
ine the convergence. 

Figure [1] shows the energy gain per site obtained by 
the AF-RVB stripe states as a function of doping den- 
sity in the t — t' — t" — J model with renormalized t*. 
The reference ground state is the uniform d-RVB wave 
function. Comparing to our previous work (l7j , the hop- 
ping modulation due to electron-phonon coupling indeed 
stabilizes the AF-RVB stripe states in the underdoped 
regime (6 ~ 0.05 - 0.15). The AF-RVB stripes with 
different periods of magnetic modulation like 8ao, 10ao, 
and 12ao have its minimum at doping density 1/8, 1/10 
and 1/12 respectively, which is exactly half-doped with 
1/2 hole per charge domain. We have not found a sta- 
ble stripe with period 6ao- Above doping density greater 
than 1/8, none of the stripe states are particularly sta- 
ble. These results are quite consistent with the "Yamada 
plot" observed in several experiments 0, 0, El, except 
that we have not studied the diagonal stripes below hole 
density of 0.05. Furthermore, to make sure 8ao is indeed 
most stable at 1/8 doping, we also calculate the energy 
of a "one-doped stripe" with spin modulation 16ao in a 
16 x 16 lattice. The AF-RVB stripe with 16ao is found 
to be energetically unstable with respect to 8ao- 

At 1/8 doping, Figure[2ja) shows that in contrast with 
6oo and 12ao, the energy gain of AF-RVB stripe state 
with 8ao increases with the strength of A. As shown in 
FiglUb), the difference in the average hole density be- 
tween hole-rich and hole-poor domains becomes larger 
as A increases. Among all possible magnetic periods of 
stripes, only the stripe with 8a,Q is stable at 1/8 dop- 
ing. We also find that long-range pair-pair correlation 
function for the AF-RVB stripe states is reduced at 1/8 
doping (not shown), which could explain the supercon- 
ducting transition temperature is suppressed near 1/8 
doping for LBCO 

In order to understand the stability of half-doped 
stripe, we examine the spatial distribution of hole den- 
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FIG. 3: (Color online) The average hole density of the AF- 
RVB stripes with the (a) 8ao, (b) 10ao (squares) and 12<zo 
(triangles) vs the doping density S using the same parame- 
ters in Fig[T] Different symbols denote different average hole 
densities along y direction. Blue dashed line is the doping 
density. 



sity rih(i)(= 1 — c^c^) as a function of doping density 
S. In FigEDJa), the 8ao stripes only has two distinct sites 
along y direction. When the difference of hole density 
between these two sites is getting larger, the state gains 
more kinetic energy from the modulated hopping terms 
as there are more holes sitting at sites with larger t*. 
Interestingly, the maximum energy gain (see the filled 
circles in FigfT])) with the largest difference in hole den- 
sity exactly happens at 1/8 doping for 8ao stripe. At 
doping density 1/6 or higher, the kinetic energy for the 
uniform state is already quite substantial, hence the sep- 
aration of hole-rich sites and hole-poor sites does not gain 
enough energy and we do not find a stable periodic stripe 
with 8ao. This argument also holds for stripe with other 
periods so there is no stable periodic stripe state in op- 
timum or overdoped systems. In the very underdoped 
regime, since the hole density is already small, separa- 
tion into a hole-rich and a hole-poor site does not bring 
a large difference to gain enough kinetic energy to stabi- 
lize the periodic stripes. However, for stripes with 10ao 
and 12ao shown in FigJ3Jb), there are three different sites 
with only one of them having more hole density than the 
uniform state. There is again the maximum difference in 
hole density occurs at smaller dopings (1/10 and 1/12, 
respectively). Similar argument also holds for lOao and 
12ao stripes. 

The argument provided above for the stability of half- 
doped stripes could be applicable to other kinds of 
stripes. Hence we also examine the site-centered AF- 
RVB stripes and antiphase stripes at 1/8 doping. As 
expected, the site-centered stripes have similar results as 
bond-centered stripes discussed above except with higher 
per-site energy (0.0022t) at 8ao- Similarly for the an- 
tiphase stripes 8ao has lowest energy at 1/8 doping but 
it again has a slightly higher per-site energy (0.0008t) 
& 

Finally, we investigate the stability of the AF-RVB 
stripes h\t — t' — t" — J system with long-range Coulomb 
interaction. In this case, the calculations omit the 
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FIG. 4: (Color online) (a) The optimized energy difference 
between the AF-RVB stripe and uniform d-RVB states as a 
function of the strength of Coulomb interaction Vc at 1/8 
doping, (b) The optimized energy difference vs doping density 
S for Vc/t = 1 



mass renormalization in the Hamiltonian for the mo- 
ment. The Coulomb interaction between holes is given 
by VcJ2i<j n h{i)nh(j)/\?i - tj\. We have only studied 
stripes with periods 6ao, 8ao and 12ao- The modulation 
of hole density and staggered magnetization is enhanced 
by Coulomb interaction. In Fig[J|a), as Vc increases, 
only the 8ao stripe pattern becomes more stable which is 
consistent with previous DMRG studies As shown 
in Fig[4|b), the 6ao and 12ao stripes always have higher 
energy than uniform d-RVB state even for Vc — t. Al- 
though the stripe with 8ao is stabilized at 1/8 doping, its 
energy gain is about half of that obtained by the mass 
renormalization effect shown in FigJ^a). These results 
indicate that although the long-range Coulomb interac- 
tion could help the formation of stripes, it is not partic- 
ularly effective in getting the half-doped stripes as the 
electron-phonon coupling. 



To summarize, we have successfully included the effect 
of mass renormalization due to weak electron-phonon in- 
teraction with the extended t — J Hamiltonian in a self- 
consistent VMC calculation. It is shown that the vertical 
half-doped stripes are most preferred in the underdoped 
regime with hole density 1/12 < <5 < 1/8. This explains 
the " Yamada plot" mostly seen in LSCO materials. Since 
the stability of stripes is enhanced by electron-phonon 
coupling, our result also implies that this coupling is 
stronger in LSCO materials than other cuprates which is 
consistent with the analysis of kink in Ref . [26| . Addition- 
ally, we have found that long-range Coulomb interaction 
does not play as an important role in the formation of 
half-doped stripes as the mass renormalization effect due 
to electron-phonon interaction. 

This work is supported by the National Science Council 
in Taiwan with Grant no.98-2112-M-001-017-MY3. The 
calculations are performed in the National Center for 
High-performance Computing in Taiwan. 
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